Chapter Five 


Controllability and Observability 


Controllability and observability represent two major concepts of modern control 
system theory. These originally theoretical concepts, introduced by R. Kalman 
in 1960, arc particularly important for practical implementations. They can be 
roughly defined as follows. 

Controllability: In order to be able to do whatever we want with the given 
dynamic system under control input, the system must be controllable. 

Observability: In order to see what is going on inside the system under 
obsen’ation, the system must be observable. 

Even though the concepts of controllability and observability are almost ab¬ 
stractly defined, we now intuitively understand their meanings. The remaining 
problem is to produce some mathematical check up tests and to define controlla¬ 
bility and observability more rigorously. Our intention is to reduce mathematical 
derivations and the number of definitions, but at the same time to derive and 
define very clearly both of them. In that respect, in Section 5.1, we start with 
observability derivations for linear discrete-time invariant systems and give the 
corresponding definition. The observability of linear discrete systems is very 
naturally introduced using only elementary linear algebra. This approach will 
be extended to continuous-time system observability, where the derivatives of 
measurements (observations) have to be used. Section 5.2. Next, in Sections 5.3 
and 5.4 we define controllability for both discrete- and continuous-time linear 
systems. 

In this chapter we show that the concepts of controllability and observability 
are related to linear systems of algebraic equations. It is well known that a 
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solvable system of linear algebraic equations has a solution if and only if the rank 
of the system matrix is full (see Appendix C). Observability and controllability 
tests will be connected to the rank tests of ceratin matrices, known as the 
controllability and observability matrices. 

At the end of this chapter, in Section 5.5, we will introduce the concepts 
of system stabilizability (detectability), which stand for controllability (observ¬ 
ability) of unstable system modes. Also, we show that both controllability and 
observability are invariant under nonsingular transformations. In addition, in the 
same section the concepts of controllability and observability are clarified using 
different canonical forms, where they become more obvious. 

The study of observability is closely related to observer (estimator) design, 
a simple, but extremely important technique used to construct another dynamic 
system, the observer (estimator), which produces estimates of the system state 
variables using information about the system inputs and outputs. The estimator 
design is presented in Section 5.6. Techniques for constructing both full-order and 
reduced-order estimators are considered. A corresponding problem to observer 
design is the so-called pole placement problem. It can be shown that for a 
controllable linear - system, the system poles (eigenvalues) can be arbitrarily 
located in the complex plane. Since this technique can be used for system linear 
feedback stabilization and for controller design purposes, it will be independently 
presented in Section 8.2. 

Several examples are included in order to demonstrate procedures for ex¬ 
amining system controllability and observability. All of them can be checked 
by MATLAB. Finally, we have designed the corresponding laboratory experi¬ 
ment by using the MATLAB package, which can contribute to better and deeper 
understanding of these important modern control concepts. 

Chapter Objectives 

This chapter introduces definitions of system controllability and observabil¬ 
ity. Testing controllability and observability is replaced by linear algebra prob¬ 
lems of finding ranks of certain matrices known as the controllability and ob¬ 
servability matrices. After mastering the above concepts and tests, students will 
be able to determine system initial conditions from system output measurements, 
under the assumption that the given system is observable. As the highlight of 
this chapter, students will learn how to construct a system's observer (estimator), 
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which for an observable system produces the estimates of state variables at any 
time instant. 

5.1 Observability of Discrete Systems 

Consider a linear, time invariant, discrete-time system in the state space form 

x(& + 1) = A d.x(k), x(0) = x 0 = unknown (5.1) 

with output measurements 

y(k) = C..x(A') (5.2) 

where x(7r) E If". y (k) E 3? p . A,/ and C d arc constant matrices of appropriate 
dimensions. The natural question to be asked is: can we learn everything about 
the dynamical behavior of the state space variables defined in (5.1) by using 
only information from the output measurements (5.2). If we know x 0 , then the 
recursion (5.1) apparently gives us complete knowledge about the state variables 
at any discrete-time instant. Thus, the only thing that we have to determine from 
the state measurements is the initial state vector x(0) = x 0 . 

Since the «-dimensional vector x(0) has n unknown components, it 
is expected that n measurements are sufficient to determine x 0 . Take 
k = 0, I ,..., v — 1 in (5.1) and (5.2), i.e. generate the following sequence 

y(0) = C d x(0) 

y(l) = C d x(l) = C d A d x(0) 

y(2) = c d x( 2 ) = C d A d x(l) = C d A^x(0) (5.3) 


y(n - 1) = C d x(n - 1) = C d A” a x(0) 


or, in matrix form 


to P 0 

(np)xl 

C d 

C d A d 

CdAj 

y( n -!). 


O 

Sc 

> •• 

SC 5 

1 


X x(0) 


(5.4) 
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We know from elementary linear algebra that the system of linear algebraic 
equations with n unknowns, (5.4), has a unique solution if and only if the system 
matrix has rank n. In this case we need 


rank 


C d 

C d A d 

C d Aj 


= n 


(5.5) 


C d A 


ri¬ 

el 


1 


Thus, the initial condition xo is completely determined if the so-called observ¬ 
ability matrix, defined by 


0{A d , C d ) 


C d 

C d A d 

C d A* 


(np)xn 


(5.6) 


C d A 


n— 

d 


1 


has rank n, that is 


rarikO = n 


(5.7) 


The previous derivations can be summarized in the following theorem. 

Theorem 5.1 The linear discrete-time system (5.1) with measurements (5.2) 
is observable if and only if the observability matrix (5.6) has rank equal to n. 

A simple second-order example demonstrates the procedure for examining 
the observability of linear discrete-time systems. More complex examples cor¬ 
responding to real physical control systems will be considered in Sections 5.7 
and 5.8. 

Example 5.1: Consider the following system with measurements 


%i(k + 1) 


T 2' 

xi (k) 

x- 2 {k + 1)_ 


3 4_ 

42 {k)\ 


y(k) = [1 


2 


x 1 (k) 

x 2 (k) 
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The observability matrix for this second-order system is given by 

1 2 ' 

7 10 _ 

Since the rows of the matrix O are linearly independent, then rank O = 2 = n, 
i.e. the system under consideration is observable. Another way to test the 
completeness of the rank of square matrices is to find their determinants. In 
this case 


O = 


C d 

C d A d 


detd? = -4 / 0 o- full rank = n = 2 


o 

Example 5.2: Consider a case of an unobservable system, which can be 


obtained by slightly modifying Example 
measurement matrices are given by 



The observability matrix is 



so that ranker = 1 < 2, and the system 


5.1. The corresponding system and 

C d = [l 2] 

2 

- IT) _ 

is unobservable. 


o 


5.2 Observability of Continuous Systems 

A linear, time invariant, continuous system in the state space form was studied 
in Chapter 3. For the puipose of studying its observability, we consider an 
input-free system 


x(t) = Ax(t), x(fo) = x 0 = unknown (5.8) 

with the corresponding measurements 


y(0 = Cx(t) 


( 5 . 9 ) 
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of dimensions x(f) G y (t) G A G 3?™ x ™, and C G 3? pX ™. Following the 
same arguments as in the previous section, we can conclude that the knowledge 
of x 0 is sufficient to determine x(f) at any time instant, since from (5.8) we have 

x(i) = e A (*-*°) x (/ 0 ) (5.10) 

The problem that we are faced with is to find x(ffi) from the available mea¬ 
surements (5.9). In Section 5.1 we have solved this problem for discrete-time 
systems by generating the sequence of measurements at discrete-time instants 
k — 0,1,2— 1, i.e. by producing relations given in (5.3). Note that a 
time shift in the discrete-time corresponds to a derivative in the continuous-time. 
Thus, an analogous technique in the continuous-time domain is obtained by taking 
derivatives of the continuous-time measurements (5.9) 

y (to) - Cx(ffi) 
y {t 0 ) - Cx(fo) = CAx(fo) 

y(/ 0 ) - Cx(ffi) - CA 2 x(/ 0 ) ( 5 . 11 ) 


y (vi ' ] (lu) = C'x f;i IJ (/.u) = CA 'x(fo) 

Our goal is to generate n linearly independent algebraic equations in n unknowns 
of the state vector x(ffi). Equations (5.11) comprise a system of np linear 
algebraic equations. They can be put in matrix form as 


yfo) 

(np)xl 

C ■ 

(np)xn 

y(fo) 


CA 


yfo) 

— 

CA 2 

x x(f 0 ) = C’x(fo) = Y (to) (5.12) 

_y (n_1) ( h). 


.CA™ -1 . 



where O is the observability matrix already defined in (5.6) and where the 
definition of Y(ffi) is obvious. Thus, the initial condition x(/ tJ ) can be determined 
uniquely from (5.12) if and only if the observability matrix has full rank, i.e. 
rankO = n. 

As expected, we have obtained the same observability result for both 
continuous- and discrete-time systems. The continuous-time observability the¬ 
orem, dual to Theorem 5.1, can be formulated as follows. 



CONTROLLABILITY AND OBSERVABILITY 


227 


Theorem 5.2 The linear continuous-time system (5.8) with measurements 
(5.9) is observable if and only if the observability matrix has full rank. 

It is important to notice that adding higher-order derivatives in (5.12) cannot 
increase the rank of the observability matrix since by the Cayley-Hamilton 
theorem (see Appendix C) for k > n we have 

71 — 1 

= (5.13) 

8 = 0 

so that the additional equations would be linearly dependent on the previously 
defined n equations (5.12). The same applies to the discrete-time domain and 
the corresponding equations given in (5.4). 

There is no need to produce a test example for the observability study of 
continuous-time systems since the procedure is basically the same as in the case 
of discrete-time systems studied in the previous section. Thus, Examples 5.1 
and 5.2 demonstrate the presented procedure in this case also; however, we have 
to keep in mind that the corresponding matrices A and C describe systems 
which operate in different time domains. Fortunately, the algebraic procedures 
are exactly the same in both cases. 


5.3 Controllability of Discrete Systems 

Consider a linear - discrete-time invariant control system defined by 

x(& + 1) = A d x(k) + Bdu(k), x(0) = x 0 (5.14) 

The system controllability is roughly defined as an ability to do whatever we 
want with our system, or in more technical terms, the ability to transfer our 
system from any initial state x(0) = x 0 to any desired final state x(£q) = xy 
in a finite time, i.e. for k\ < oo (it makes no sense to achieve that goal at 
k\ = oo). Thus, the question to be answered is: can we find a control sequence 
u(0), u(l),..., u(n — 1), such that x(fc) = xy? 

Let us start with a simplified problem, namely let us assume that the input 
u (k) is a scalar, i.e. the input matrix B ( / is a vector denoted by hj. Thus, we have 


x(k + 1) = A dx(k) + bd'u(k), x(0) = x, 


(5.15) 
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Taking k — 0,1,2,.... ri in (5.15), we obtain the following set of equations 
x(l) = Arfx(O) + b d 'tt(0) 

x(2) = A d x(l) + b d u(l) - A^x(0) + A d b d w(0) + b d u(l) 

(5.16) 

x( ri) - A^x(0) + A^ _1 b c i'u(0) +-f- b d u(n - 1) 

The last equation in (5.16) can be written in matrix form as 

u{n — 1) 
u(n - 2) 

x(»t) - A;;x(0) - b d : A d b d : • • • : A^ _1 b d : (5.17) 

L ' «(!) 

_ tt(0) . 

Note that bdiA^bd: • • • :A) _ 1 b,j is a square matrix. We call it the controlla¬ 
bility matrix and denote it by C. If the controllability matrix C is nonsingular, 
equation (5.17) produces the unique solution for the input sequence given by 

u(n — 1) 
u(n — 2) 

: — C 1 (x(/q) — A;.x(()).| (5.18) 

“(1) 

_ «(0) . 

Thus, for any x(r«) = xy, the expression (5.18) determines the input sequence 
that transfers the initial state x 0 to the desired state xy in n steps. It follows 
that the controllability condition, in this case, is equivalent to nonsingularity of 
the controllability matrix C. 

In a general case, when the input u (k) is a vector of dimension r, the 
repetition of the same procedure as in (5.15)—(5.17) leads to 

u (n — 1) 
u(n - 2) 

x(n)-A3x(0)= BjiAdBdi iApBd | (5.19) 

u(l) 

_ u(0) _ 
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The controllability matrix, in the general vector input case, defined by 


C(A d ,B d ) 


B d :A d B d : ••• : 


(5.20) 


is of dimension n X r ■ n. The corresponding system of n linear alge¬ 
braic equations in r ■ n unknowns for n r -dimensional vector components of 
u(0), n( I). u(n - 1), given by 


Qnx(nr) 


u (n — 1) 
u(n - 2) 


(nr)xl 


u (l) 

u(0) 


x(n) - A.:x(0) 


x / - A d x (0) 


(5.21) 


will have a solution for any x/ if and only if the matrix C has full rank, i.e. 
rankC = n (see Appendix C). 

The controllability theorem is as follows. 

Theorem 5.3 The linear discrete-time system (5.14) is controllable if and 
only if 


rankC = n (5.22) 

where the controllability matrix C is defined by (5.20). 


5.4 Controllability of Continuous Systems 

Studying the concept of controllability in the continuous-time domain is more 
challenging than in the discrete-time domain. At the beginning of this section we 
will first apply the same strategy as in Section 5.3 in order to indicate difficulties 
that we arc faced with in the continuous-time domain. Then, we will show how 
to find a control input that will transfer our system from any initial state to any 
final state. 

A lineal - continuous-time system with a scalar input is represented by 

x = Ax+ b«, x(t 0 ) = x 0 (5.23) 
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Following the discussion and derivations from Section 5.3, we have, for a scalar 
input, the following set of equations 

• d 

X = —X = Ax + DU 

at 

d 2 

x = ——x = A 2 x + Abu + bit 

dt 2 T T (5.24) 


d n 

X W = x = A”x + A^'bu + A n ~ 2 hu + • • • + biM-V 
dt n 

The last equation in (5.24) can be written as 

u^ n ~ 2 \t) 

xW(i)-A n x(t) = C • (5.25) 

u(t) 

. u (4 . 

Note that (5.25) is valid for any t 6 [ta,tf] with // free but finite. Thus, the 
nonsingularity of the controllability matrix C implies the existence of the scalar 
input function u( l) and its n — 1 derivatives, for any t < tf < oo. 

For a vector input system dual to (5.23), the above discussion produces the 
same relation as (5.25) with the controllability matrix C given by (5.20) and with 
the input vector u(/ ) e f?'. that is 

u( n - 2 )(/) 

C nxm ’ n j = xW(i)-A n x(t) = 7 (t) (5.26) 

“(*) 

. u (0 . 

It is well known from linear algebra that in order to have a solution of (5.26), 
it is sufficient that 


rankC = rank C : 7 (t) 


(5.27) 
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Also, a solution of (5.26) exists for any 7 (f)—any desired state at t —if and only if 

rankC = n (5.28) 


Equations (5.25) and (5.26) establish relationships between the state and 
control variables. However, from (5.25) and (5.26) we do not have an explicit 
answer about a control function that is transferring the system from any initial 
state x((q) to any final state x(t 1 ) = xj. Thus, elegant and simple derivations 
for the discrete-time controllability problem cannot be completely extended to 
the continuous-time domain. Another approach, which is mathematically more 
complex, is required in this case. It will be presented in the remaining part of 
this section. 

From Section 3.2 we know that the state space equation with the control 
input has the following solution 

t 

x(t) = e A ^~ io ^x(to) + J e A ( f-T )Bu(r)dr 
*0 

At the final time t\ we have 

t\ 

x(ti ) = xj = e A ^ 1 _t °'x(fo) + J e A ( tl_T )Bu(r)dr 

to 


e~ At 'x f - e- Aio x(t 0 ) = J e“ AT Bu(r)dr 

to 

Using the Cayley-Hamilton theorem (see Appendix C), that is 

n-l 

e- AT = ^a ! -(r)A' (5.29) 

i=0 

where a^(r), i = 0 , 1 , — 1 , are scalar time functions, the previous equation 

can be rewritten as 


n — 1 

e~ Atx xj - e~ At °x{t 0 ) = A B 

z=0 


J «i(r)u(r )dr 

to 
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or 


e ‘x 


/ - e- At °x(*o) = 


B:AB: ••• iA n_1 B 


r 1 1 

/ a 0 ('r)u(r)dr 

to 
t j. 

J ai(r)u(r)dr 
^0 


/a n _i(r)u(r)dr 

L 0 


On the left-hand side of this equation all quantities are known, i.e. we have a 
constant vector. On the right-hand side the controllability matrix is multiplied 
by a vector whose components are functions of the required control input. Thus, 
we have a functional equation in the form 


const”* 1 =C( A,B)” X ™ 


fiKr)) 

f 2 (u(r)) 


t e 


Lf n _i(u(r)) J 


(5.30) 


A solution of this equation exists if and only if rankC(A,B) = re, which is the 
condition already established in (5.28). In general, it is very hard to solve this 
equation. One of the many possible solutions of (5.30) will be given in Section 
5.8 in terms of the controllability Grammian. The controllability Grammian is 
defined by the following integral 


^1 

W(t 0 ,*i) = J e A( ‘ 0_T) BB T e AT(t °- T UT 

to 


(5.31) 


The results presented in this section can be summarized in the following 
theorem. 

Theorem 5.4 The linear continuous-time system is controllable if and only if 
the controllability matrix C has full rank, i.e. rankC = n. 

We have seen that controllability of linear continuous- and discrete-time 
systems is given in terms of the controllability matrix (5.20). Examining the 
rank of the controllability matrix comprises an algebraic criterion for testing 
system controllability. The example below demonstrates this procedure. 
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Example 5.3: Given the linear continuous-time system 


X = 

' 0 

3 

1 

-4 

-2' 

5 

x + 

'0 

2 

-1' 

-3 


-6 

7 

8 


4 

-5 


The controllability matrix for this third-order system is given by 


C 


b ; ab :a 2 b 


0 

-1 

-6 

7 


2 

-3 

12 

-10 

A 2 B 

4 

-5 

46 

-55 



Since the first three columns arc linearly independent we can conclude that 
rankC = 3. Hence there is no need to compute A 2 B since it is well known 
from lineal - algebra that the row rank of the given matrix is equal to its column 
rank. Thus, rankC = 3 = n implies that the system under consideration is 
controllable. 


o 


5.5 Additional Controllability/Observability Topics 

In this section we will present several interesting and important results related to 
system controllability and observability. 

Invariance Under Nonsingular Transformations 

In Section 3.4 we introduced the similarity transformation that transforms 
a given system from one set of coordinates to another. Now we will show 
that both system controllability and observability are invariant under similarity 
transformation. 

Consider the vector input form of (5.23) and the similarity transformation 

x = Px (5.32) 


x = Ax + Bu 


such that 
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where A = PAP 1 and B = PB. Then the following theorem holds. 

Theorem 5.5 The pair (A, B) is controllable if and only if the pair ^A, B 
is controllable. 

This theorem can be proved as follows 


BAB: 

• A^'B 


PB:PAP 

X PB : • • • 

:PA" _1 P _1 PB 

P B:AB: 

• • : A 7i_1 B = PC( A, B) 


Since P is a nonsingular matrix (it cannot change the rank of the product PC), 
we get 


rankC^A,B^ =rankC(A,B) 


which proves the theorem and establishes controllability invariance under a 
similarity transformation. 

A similar theorem is valid for observability. The similarity transformation 
(5.32) applied to (5.8) and (5.9) produces 

x = Ax 
y = Cx 

where 

C = CP 1 


Then, we have the following theorem 

Theorem 5.6 The pair (A, C) is observable if and only if the pair ^A, C 
is observable. 


The proof of this theorem is as follows 


c(a,c) = 

■ C ■ 

CA 

CA 2 

— 

CP 1 

CP X PAP 1 
CP X PA 2 P 1 

— 

c 

CA 

CA 2 


.CA™- 1 . 


.CP^PA™ X P 1 . 


.CA™- 1 . 
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that is, 

e>(A,c) = 0(A,C)P _1 

The nonsingularity of P implies 

rankC7^A,C^ =rankC?(A,C) 

which proves the stated observability invariance. 

Note that Theorems 5.5 and 5.6 are applicable to both continuous- and 
discrete-time linear systems. 

Frequency Domain Controllability and Observability Test 

Controllability and observability have been introduced in the state space 
domain as pure time domain concepts. It is interesting to point out that in the 
frequency domain there exists a very powerful and simple theorem that gives a 
single condition for both the controllability and the observability of a system. 
It is given below. 

Let H(s) be the transfer function of a single-input single-output system 

H(s ) = c(al - A) _1 b 

Note that H(s) is defined by a ratio of two polynomials containing the cor¬ 
responding system poles and zeros. The following controllability-observability 
theorem is given without a proof. 

Theorem 5.7 If there are no zero-pole cancellations in the transfer function 
of a single-input single-output system, then the system is both controllable and 
observable. If the zero-pole cancellation occurs in H (s), then the system is either 
uncontrollable or unobservable or both uncontrollable and unobservable. 

A similar theorem can be formulated for discrete linear time invariant 
systems. 

Example 5.4: Consider a linear continuous-time dynamic system represented 
by its transfer function 

(s + 3) _ s + 3 

(7 + l)(s + 2)(s + 3) ~ s 3 + 6s 2 + 11s + 6 

Theorem 5.7 indicates that any state space model for this system is either 
uncontrollable or/and unobservable. To get the complete answer we have to go 
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to a state space form and examine the controllability and observability matrices. 
One of the possible many state space forms of II (s) is as follows 


Xi ' 


'-6 -11 -6' 


~X\ 

X2 

= 

1 0 0 


X-2 

X3_ 


0 1 0 _ 


X-.i 



X 

ll 


y = [ o i 3] 

X2 



. X 3. 


It is easy to show that the controllability and observability matrices arc given by 


Since 


'1 

-6 

25 ' 


' 0 

1 

3 ' 

0 

1 

-6 

, o = 

1 

3 

0 

_0 

0 

1 


-3 

-11 

-6 


det C = 1 j: 0 => rankC = 3 = n 


and 


det 0 = 0 rankO <3 = n 


this system is controllable, but unobservable. 

Note that, due to a zero-pole cancellation at s = —3, the system transfer 
function H{s) is reducible to 


H(s) = H r (s ) 


1 

(s + l)(s + 2) 


1 

s 2 + 3s + 2 


so that the equivalent system of order n 
form _ 


i 


1 

to 

-31 

[X2 r. 


Li 

-1 

O 


V — [0 1 ] 

For this reduced-order system we have 


2 has the coiTesponding state space 


X\ r 


'i' 


+ 

o_ 

U 

_ x 2r . 




Xir 

p2r 


c 


1 -2 
0 1 ’ 


o 


0 1 
1 0 
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and therefore the system is both controllable and observable. 

Interestingly enough, the last two mathematical models of dynamic systems 
of order n — 3 and n — 2 represent exactly the same physical system. Apparently, 
the second one (n — 2) is preferred since it can be realized with only two 
integrators. 


o 

It can be concluded from Example 5.4 that Theorem 5.7 gives an answer 
to the problem of dynamic system reducibility. It follows that a single-input 
single-output dynamic system is irreducible if and only if it is both controllable 
and observable. Such a system realization is called the minimal realization. If 
the system is either uncontrollable and/or unobservable it can be represented 
by a system whose order has been reduced by removing uncontrollable and/or 
unobservable modes. It can be seen from Example 5.4 that the reduced system 
with n — 2 is both controllable and observable, and hence it cannot be further 
reduced. This is also obvious from the transfer function H r (s). 

Theorem 5.7 can be generalized to multi-input multi-output systems, where 
it plays very important role in the procedure of testing whether or not a given 
system is in the minimal realization form. The procedure requires the notion of 
the characteristic polynomial for proper rational matrices which is beyond the 
scope of this book. Interested readers may find all details and definitions in 
Chen (1984). 

It is important to point out that the similarity transformation does not change 
the transfer function as was shown in Section 3.4. 

Controllability and Observability of Special Forms 

In some cases, it is easy to draw conclusions about system controllability 
and/or observability by examining directly the state space equations. In those 
cases there is no need to find the corresponding controllability and observability 
matrices and check their ranks. 

Consider the phase variable canonical form with 

x = Ax + Bu 
y = Cx 

where 
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' 0 

1 

0 

0 ' 


TT 

A = 

0 

0 

1 

0 

, B = 

0 


0 

0 

0 

1 


0 


-a 0 

-fll 

-a 2 ■ ■ 

1 

© 

1 


.1. 


C = [1 0 0 ... 0] 

This form is both controllable and observable due to an elegant chain connection 
of the state variables. The variable #i (t) is directly measured, so that x 2 {t) 
is known from x 2 {t) = Also, x^{t) = x 2 {t) = x\{t), and so on, 

x n (t ) = x ^(f). Thus, this form is observable. The controllability follows 
from the fact that all state variables arc affected by the control input, i.e. x n is 
affected directly by u{t) and then x n -\(t) by x n (u(t)) and so on. The control 
input is able to indirectly move all state variables into the desired positions so 
that the system is controllable. This can be formally verified by forming the 
corresponding controllability matrix and checking its rank. This is left as an 
exercise for students (see Problem 5.13). 

Another example is the modal canonical form. Assuming that all eigenvalues 
of the system matrix arc distinct, we have 

x = AxT Tu 

y = 3x 

where 



-Ar 

0 • 

• • 0 - 


'71' 

A = 

0 

a 2 • 

• • 0 

, r = 

72 


.0 

0 • 

^ n - 


- In - 


3 = [Si S 2 • • • 8 n ] 

We arc apparently faced with n completely decoupled first-order systems. Ob¬ 
viously, for controllability all 7 i, i = 1 ...., n. must be different from zero, 
so that each state variable can be controlled by the input u(f). Similarly, 
/ 0 , i = 1 ,..., n, ensures observability since, due to the state decomposition, 
each system must be observed independently. 
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The Role of Observability in Analog Computer Simulation 

In addition to applications in control system theory and practice, the concept 
of observability is useful for analog computer simulation. Consider the problem 
of solving an nth-order differential equation given by 

n m 

y {n) + E a n-iy {h ~ l) = E 

i—1 z=0 

with known initial conditions for .(/(()). /)(()),.... **(()). This system can be 

solved by an analog computer by using n integrators. The outputs of these n 
integrators represent the state variables xi, x %,..., x n , so that this system has the 
state space form 


x = Ax + bn, x(0) = unknown 

y = cx 

However, the initial condition for x(0) is not given. In other words, the 
initial conditions for the considered system of n integrators arc unknown. They 
can be determined from j/(0), r/(0),..., i/( n_1 )(0) by following the observability 
derivations performed in Section 5.2, namely 

2/(0) = cx(0) 

.(/(()) = cx(0) = cAx(O) + cbn(0) 

j/(0) = cx(0) = cA 2 x(0) + cAbn(O) + cbn(0) 

.(/('^'>(()) = cx^-O(O) = cA n_1 x(0) + cA Tl_2 bu(0) 

+ cA Tl_3 bn(0) + • • • + cAbi/ 1 ” 3i (0) + cbn f:<i ' 2 >(()) 

This system can be written in matrix form as follows 


1 

cT o' . . 

i_ 

= O ■ x(0) + T 

1 

£ 

" 3 ° 

_ I 

i 

i 

O 

1_ 


.«("- _2) (0). 


where O is the observability matrix and T is a known matrix. Since 
if(0), «(0),..., u(" -1 )(0) arc known, it follows that a unique solution for x(0) 



240 


CONTROLLABILITY AND OBSERVABILITY 


exists if and only if the observability matrix, which is square in this case, is 
invertible, i.e. the pair (A,c) is observable. 

Example 5.5: Consider a system represented by the differential equation 


W + 4 Jt +4y= Ji + U: = 2; -^ () ) '• L u i t ) = e 4t , l >° 


Its state space form is given by 

x = Ax + bu = 


0 1 
-4 -4 


x 


y = cx = [1 1 ]x 

The initial condition for the state space variables is obtained from (5.33) as 


leading to 



y(o) 

j(o) 


o 

cbu(O) 


2 

1 


0 

1 


1 

-4 



2 

0 


x(0) 


>i(0)' 


'-6' 

*2(0). 


. 8 . 


This means that if analog computer simulation is used to solve the above second- 
order differential equation, the initial conditions for integrators should be set to 
-6 and 8. 


o 


Stabilizability and Detectability 

So far we have defined and studied observability and controllability of the 
complete state vector. We have seen that the system is controllable (observable) 
if all components of the state vector are controllable (observable). The natural 
question to be asked is: do we really need to control and observe all state 
variables? In some applications, it is sufficient to take care only of the unstable 
components of the state vector. This leads to the definition of stabilizability and 
detectability. 

Definition 5.1 A linear system (continuous or discrete) is stabilizable if ail 
unstable modes are controllable. 
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Definition 5.2 A linear system (continuous or discrete) is detectable if all 
unstable modes are obsen’able. 

The concepts of stabilizability and detectability play very important roles in 
optimal control theory, and hence are studied in detail in advanced control theory 
courses. For the puipose of this course, it is enough to know their meanings. 

5.6 Observer (Estimator) Design 1 

Sometimes all state space variables are not available for measurements, or it is not 
practical to measure all of them, or it is too expensive to measure all state space 
variables. In order to be able to apply the state feedback control to a system, 
all of its state space variables must be available at all times. Also, in some 
control system applications, one is interested in having information about system 
state space variables at any time instant. Thus, one is faced with the problem 
of estimating system state space variables. This can be done by constructing 
another dynamical system called the observer or estimator, connected to the 
system under consideration, whose role is to produce good estimates of the state 
space variables of the original system. 

The theory of observers stalled with the work of Luenberger (1964, 1966, 
1971) so that observers are very often called Luenberger observers. According 
to Luenberger, any system driven by the output of the given system can serve 
as an observer for that system. Two main techniques are available for observer 
design. The first one is used for the full-order observer design and produces 
an observer that has the same dimension as the original system. The second 
technique exploits the knowledge of some state space variables available through 
the output algebraic equation (system measurements) so that a reduced-order 
dynamic system (observer) is constructed only for estimating state space variables 
that are not directly obtainable from the system measurements. 

5.6.1 Full-Order Observer Design 

Consider a linear time invariant continuous system 

x(t) = Axffi + Bu(f), x(/n) = x n = unknown 

(5 34) 

y (t) = Cx(t) 1 ; 

1 This section may be skipped without loss of continuity. 
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where x e 9?”, u E 3?’’, y e with constant matrices A,B,C having 
appropriate dimensions. Since from the system (5.34) only the output variables, 
y (t), arc available at all times, we may construct another artificial dynamic 
system of order n (built, for example, of capacitors and resistors) having the 
same matrices A, B, C 


x(t) - Ax(t) + Bu(f), x(t 0 )-x o 

y(0 = c ±(t) 


(5.35) 


and compare the outputs y(/) and y{i). Of course these two outputs will be 
different since in the first case the system's initial condition is unknown, and in 
the second case it has been chosen arbitrarily. The difference between these two 
outputs will generate an error signal 


y (t) - y(t ) = C x(t) - Cx(t) - Ce(t) (5.36) 


which can be used as the feedback signal to the artificial system such that the 
estimation (observation) error e(t ) = x(t) - x(t) is reduced as much as possible. 
This can be physically realized by proposing the system-observer structure as 
given in Figure 5.1. 



Figure 5.1: System-observer structure 
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In this structure K represents the observer gain and has to be chosen such 
that the observation error is minimized. The observer alone from Figure 5.1 is 
given by 

x(t) = Ax(t) + Bu(jf) + K(y(f) - y(/)) = Ax(t) + Bu({) + KCe(f) (5.37) 

From (5.34) and (5.37) it is easy to derive an expression for dynamics of the 
estimation (observation) error as 

e(f) = (A - KC)e(i) (5.38) 

If the observer gain K is chosen such that the feedback matrix A — KC is 
asymptotically stable, then the estimation error e(7) will decay to zero for any 
initial condition e(fo)» This can be achieved if the pair (A, C) is observable. 
More precisely, by taking the transpose of the estimation error feedback matrix, 
i.e. A t — C T K. T , we see that if the pair (A. T , C T ) is controllable, then we 
can do whatever we want with the system, and thus we can locate its poles in 
arbitrarily asymptotically stable positions. Note that controllability of the pair 
(A 7 . C T ) is equal to observability of the pair (A,C), see expressions for the 
observability and controllability matrices. 

In practice the observer poles should be chosen to be about ten times faster 
than the system poles. This can be achieved by setting the minimal real part of 
observer eigenvalues to be ten times bigger than the maximal real part of system 
eigenvalues, that is 

}lo&se^stev 10 

Theoretically, the observer can be made arbitrarily fast by pushing its eigenvalues 
far to the left in the complex plane, but very fast observers generate noise in 
the system. A procedure suggesting an efficient choice of the observer initial 
condition is discussed in Johnson (1988). 

It is important to point out that the system-observer structure preserves the 
closed-loop system poles that would have been obtained if the linear perfect 
state feedback control had been used. The system (5.34) under the perfect state 
feedback control, i.e. u (t) = —Fx(7) has the closed-loop form as 


x(t) = (A - BF)x(t) 


(5.39) 
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so that the eigenvalues of the matrix A — BF arc the closed-loop system poles 
under perfect state feedback. In the case of the system-observer structure, as 
given in Figure 5.1, we see that the actual control applied to both the system 
and the observer is given by 


u(f) = -Fx(i) = -F x(/) + Fe(f) (5.40) 


so that from (5.34) and (5.38) we have 


X 


A - BF 

BF 

X 

e 


0 

A - KC 

e 


(5.41) 


Since the state matrix of this system is upper block triangular, its eigenvalues are 
equal to the eigenvalues of matrices A — BF and A — KC. A very simple relation 
among x, e, and x can be written from the definition of the estimation error as 


X 


i 

0 ' 


X 

= T 

X 

e 


i 

-I 


X 


X 


(5.42) 


Note that the matrix T is nonsingular. In order to go from xe-coordinates to 
xx-coordinates we have to use the similarity transformation defined in (5.42), 
which by the main property of the similarity transformation indicates that the 
same eigenvalues, i.e. A(A — BF) and A(A — KC), are obtained in the xx- 
coordinates. 

This important observation that the system-observer configuration has closed- 
loop poles separated into the original system closed-loop poles under perfect state 
feedback and the actual observer closed-loop poles is known as the separation 
principle. 


5.6.2 Reduced-Order Observer (Estimator) 

In this section we show how to construct an observer of reduced dimension 
by exploiting knowledge of the output measurement equation. Assume that the 
output matrix C has rank l, which means that the output equation represents l 
linearly independent algebraic equations. Thus, equation 

y (t) = C x(t) (5.43) 

produces l equations for n unknowns of the state space vector x(t). Our goal 
is to construct an observer of order n — l for estimation of the remaining n — l 
state space variables. 
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The reduced-order observer design will be presented according to the results 
of Cumming (1969) and Gopinath (1968, 1971). The procedure for obtaining 
this observer is not unique, which is obvious from the next step. Assume that 
a matrix Ci exists such that 


rank 


C 

Cl 


= n 


and introduce a vector p e 9?* as 


p(/) = C 1 x(() 


(5.44) 


(5.45) 


From equations (5.43) and (5.45) we have 


x(f) 


c 

-1 

|y (ty 

Cl 


pW. 


(5.46) 


Since the vector p(f) is unknown, we will construct an observer to estimate it. 
Introduce the notation 

-l 

= [Li L 2 ] (5.47) 


C 

Cl 


so that from (5.46) we get 


x(/) = Liy(f) + L 2 pi I) 


(5.48) 


An observer for p(f) can be constructed by finding first a differential equation 
for p( t) from (5.45), that is 


p = Cix = CiAx + CiBu = CiAL 2 p + CiALiy + CiBu (5.49) 


Note that from (5.49) we arc not able to construct an observer for p(7 ) since y(7) 
does not contain explicit information about the vector p(f), but if we differentiate 
the output variable we get from (5.34) and (5.48) 


y = Cx = CAx + CBu = CAL 2 p + CALiy + CBu (5.50) 


i.e. y (t) caiiies information about p(/). 
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An observer for p(7), according to the observer structure given in (5.37), is 
obtained from the last two equations as 

p = CiAL 2 p + Ci ALiy + CiBu + K x (y - y) (5.51) 

where Ki is the observer gain. If in equation (5.50) we replace p( /) by its 
estimate, we will have 


y = CAL 2 p + CALiy + CBu (5.52) 


so that 

p = C 1 AL 2 P + CiALiy + CiBu + Ki(y — CAL 2 P — CALiy — CBu) 

(5.53) 

Since it is impractical and undesirable to differentiate y (t) in order to get y(7) 
(this operation introduces noise in practice), we take the change of variables 

q = P - Kiy (5.54) 

This leads to an observer for q(7) of the form 

q(7) = A,,q(/) + B q u(t) + K q y (t) (5.55) 


where 


A, = CiAL 2 - KiCAL 2 , B„ = CiB-KiCB 

(5.56) 

K, = CiAL 2 Ki + C 1 AL 1 - K 1 CAL 1 - KiCAL 2 Ki 

It is left as an exercise to students (see Problem 5.18) to derive (5.55) and (5.56). 
The estimates of the original system state space variables are now obtained from 
(5.48) and (5.53) as 

x(f) = Liy(2) + L 2 p(f) = L 2 q(i) + (L x + L 2 Ki)y (5.57) 

The obtained system-reduced-observer structure is presented in Figure 5.2. 
There are other ways of constructing the system observers (Luenberger, 1971; 
Chen, 1984). The reader particularly interested in observers is referred to a 
specialized book on observers for linear systems (O’Reilly, 1983). 
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Figure 5.2: System-reduced-observer structure 


5.7 MATLAB Case Study: F-8 Aircraft 


In the case of high-order systems (n > 3), obtaining the controllability and ob¬ 
servability matrices is computationally very involved. The MATLAB package for 
computer-aided control system design and its CONTROL toolbox help to over¬ 
come this problem. Moreover, use of MATLAB enables a deeper understanding 
of controllability and observability concepts. Consider the following fourth-order 
model of an F-8 aircraft studied in Teneketzis and Sandell (1977), Khalil and 
Gajic (1984), Gajic and Shen (1993). The aircraft dynamics in continuous-time 
is described by the following matrices 


-0.0135700 

-32.2 

-46.300 

0.0000 

0.0001200 

0.0 

1.214 

0.0000 

-0.0001212 

0.0 

-1.214 

1.0000 

0.0005700 

0.0 

-9.010 

-0.6696 


B = [-0.433 0.1394 -0.1394 -0.1577] r 


0 0 0 1 
10 0 0 


c = 
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By using the MATLAB function ctrb (for calculation of the controllability 
matrix C) and obsv (for calculation of the observability matrix O), it can be 
verified that this system is both controllable and observable, namely 

rankC 4 * 4 - 4, rankO 8 * 4 = 4 

By using the MATLAB function det (to calculate a matrix determinant), we get 

detC = -6.8690 

Since detC is far from zero it seems that this system is well controllable (the 
controllability margin is big). 

If we discretize the continuous-time matrices A and B by using the sampling 
period AT = 0.1, we get a somewhat surprising result. Namely 

detC(A d ,B d )| At=0 .i = -4.5381 x lO” 10 

Thus, this discrete system is almost uncontrollable. Theoretically, it is still 
controllable but we need an enormous amount of energy in order to control it. 
For example, let the initial condition be x(0) = [1 1 1 1] and let the final 

state x(4) be the coordinate origin. Then, by (5.18), the control sequence that 
solves the problem of transferring the system from x(0) to x(4) = 0 , obtained 
by using MATLAB is 


'“(!)' 


' 3.1682 ' 

u( 2) 

= 10 7 x 

-9.7464 

«( 3 ) 

9.9944 

“( 4 ). 


. 3.4163 . 


Apparently, this result is unacceptable and this discrete system is practically 
uncontrollable. 

Note that the eigenvalues of the continuous-time controllability Gram- 
mian (5.31), obtained by using the MATLAB function gram, have values 
7.48 x 10 4 , 0.42, 0.037, 0.0068. The eigenvalues of the controllability Gram- 
mian are the best indicators of the controllability measure. Since two of them 
are very close to zero, the original system is very badly conditioned from the 
controllability point of view even though detC(A,B) is far from zero. The 
interested reader can find more about controllability and observability measures 
in a very comprehensive paper by Muller and Weber (1972). 



CONTROLLABILITY AND OBSERVABILITY 


249 


5.8 MATLAB Laboratory Experiment 

Part 1. The controllability staircase form of the system 


x = Ax + Bu 
y = Cx 

clearly distinguishes controllable and uncontrollable parts of a control system. It 
can be obtained by the similarity transformation, and is defined by 


X c 


A c 

Al2 


n 

_i 

X-nc 


0 

A 

■^nc 


^nc J 


y = [C c c nc ] 



(5.58) 


where x c arc controllable modes, and x nc arc uncontrollable modes. Apparently, 
in this structure the input u cannot influence the state variables x nc ; hence these 
are uncontrollable. Similarly, one can define the observability staircase form as 


x 0 

Xy,_o 


A 0 

A 21 



x 0 

^■no 


y = [c 0 o] 


B 0 


u 


(5.59) 


with x 0 observable and x no unobservable. Due to the fact that only x 0 appeal's 
in the output and that x 0 and x no are not coupled through the state equations, 
the state variables x no cannot be observed. 

Use the MATLAB functions ctrbf (controllable staircase form) and obsvf 
(observable staircase form) to get the corresponding forms for the system 


■ l 

to 

3 

4 

5 

6 ' 


' 1 

1 ' 

-l 

-2 

-3 

-4 

-5 

-6 


-1 

1 

0 

0 

0 

0 

0 

0 

, B = 

1 

-1 

l 

1 

1 

1 

1 

1 

0 

0 

to 

2 

2 

2 

2 

2 


1 

2 

L-i 

0 

2 

-1 

0 

2 . 


_-2 

1 


C = [—1 1 -1 0 -2 2] 

Identify the corresponding similarity transformation. 
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Part 2. Derive analytically that the transfer function of (5.58) is given in 
terms of the controllable parts, i.e. it is equal to 

H c (s) = C c (sl - A C ) -1 B C = LI(>) = C(sl - A) _1 B (5.60) 

Clarify your answer by using the MATLAB function for the transfer function 
ss2zp, i.e. show that both transfer functions have the same gains, poles, and 
zeros (subject to zero-pole cancellation). 

Do the same for the observable staircase form, i.e. show that 

H 0 (s) = C 0 (sl — A 0 ) -1 B 0 = H(s) (5.61) 

and justify this identity by using the MATLAB function ss2zp. 

Part 3. Examine the controllability and observability of the power system 
composed of two interconnected areas considered in Geromel and Peres (1985) 
and Shen and Gajic (1990) 

-0 0.55 0 0 0 -5.5 0 0 0 ' 

00 1 000 0 00 

0 -3.3 -0.05 6 0 3.3 0 0 0 

0 0 0 -3.3 3.3 0 0 0 0 

A = 0 0 -5.2 0 -13 0 0 0 0 

00 0 000 1 00 

0 3.3 0 0 0 -3.3 -0.05 6 0 

0 0 0 0 0 0 0 -3.3 3.3 

.0 0 0 0 0 0 -5.2 0 -13. 

TO 000 13 000 0 1 T 

“ [o 0 0 0 0 0 0 0 13 


' 1 0.43 0 0 0 0 0 0 O' 

0 0 010 0 000 

-1 0 0 0 0 0.43 0 0 0 

.0 0 000 0 01 0 . 


Part 4. Follow the steps used in Section 5.7, but this time for the F-15 
aircraft, whose state space model was presented in Example 1.4. Consider both 
the subsonic and supersonic flight conditions. Comment on the results obtained. 
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Part 5. The controllability Grammian is defined in (5.31) as 

^1 

W(fo, ti) = J e A(fo “ r) BB T e AT(io “ r) dr 

to 

(a) Show analytically that the control input given by 

u( t ) = B 7 V A ‘ ?! W x(t u ) - ( ) (5.62) 

will drive any initial state x(f 0 ) into any desired final state x( /i) = xy. Note that 
under the controllability assumption many control inputs can be found to transfer 
the system from the initial to the final state. The expression given in (5.62) is 
also known as the minimum energy control (Klamka, 1991) since in addition to 
driving the system from x(l 0 ) to x(/ 1 ) = xj, it also minimizes an integral of the 
square of the input (energy), u T (f)u(t), in the time interval (/ U: fil¬ 
th) Using the MATLAB function gram, find the controllability Grammian 
for the system defined in Part 4 for t 0 = 0 and t\ = 1. One of several 
known controllability tests states that the system is controllable if and only if its 
controllability Grammian is positive definite (Chen, 1984; Klamka, 1991). Verify 
whether or not the controllability Grammian for this problem is positive definite. 

(c) Find the control input u(/) that drives the system defined in Paid 4 
from the initial condition x(0) = [0 0 0 0 0] 7 to the final state x(l) = 

[1111 if. 

Part 6. By duality to the controllability Grammian, the observability 
Grammian is defined as 

^1 

V(lu.h ) = J < AUiu r) C 7 C< A(?u r hlr (5.63) 

*0 

Note that the observability Grammian is in general a positive semidefinite matrix. 
It is known in the literature on observability that if and only if the observability 
Grammian is positive definite, the system is observable (Chen, 1984). Check the 
observability of the system given in Part 3 by using the observability Grammian 
test. 
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5.10 Problems 

5.1 Test the controllability and observability of the following systems 


C = [0 -1] 


'2 

-1 

0 


T—1 

1 

0 

1 

0 

, B = 

0 

1 

0 

- 1 _ 


_ 1 _ 


5.2 Find the values for parameters b\, b- 2 , and 63 such that the given system 
is controllable 


' 1 

0 

-L 


~bi 

0 

0 

2 

0 

, B = 

0 

b 2 

1 

i 

to 

0 

3 


CO 

1 

O 



5.3 Find the values for parameters c\ and c 2 such that the following system 
is observable 


-1 1 
1 -1 


A = 


C = [ci c 2 ] 
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If the output vector of the corresponding discrete system is given by 
[ 3 /( 0 ) 2/(1)] = [1 2], find the system's initial condition. 

5.4 Verify that all columns of the matrix 

' 0 1 0] 3 

A 3 =-231 
-1 0 1_ 

can be expressed as a linear combination of the columns forming matrices 
I, A, and A 2 (see 5.13). 

5.5 Assuming that the desired final state of a discrete system represented by 

0 1 °1 I" 0 ! I" 1 " 

A = -2 .3 1 , B = 1 , x(0) = 1 

-1 0 lj |_2 J |_1_ 

is x(3) = [0 —1 I ] T find the control sequence that transfers the system 

from x(0) to x(3). 

5.6 Find a solution to Problem 5.5 in the case of a two-input system that has 
the input matrix 

' 0 1 ' 

B = -1 0 

. 0 2 _ 

The remaining elements arc the same as in Problem 5.5. 

5.7 Determine conditions on b\,b 2 ,ci, and c 2 such that the following system 
is both controllable and observable 

A= E J]- B= fc]- C=[ci 

Assume that the input to this system is known. Find the initial conditions 
of this system in terms of the given input in the case when the measured 

output is y(t ) = e~ f cos(/). 

5.8 Using the frequency domain criterion, check the joint controllability and 
observability of the system 



' 1 

0 

2' 


' 1 

A = 

-1 

0 

0 

, B = 

0 


_ 1 

2 

0 


-1_ 
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5.9 Find the initial conditions of all integrators in an analog computer simulation 
of the following differential equation 

(/-’ y .. da d 2 u (l u 


d 3 


y 

dt 3 


dt 2 


3 ^ 

dt 


y = 


dt 2 dt 


2 /( 0 ) 


dyj 0) 
dt 


d 2 y( 0) 
dt 2 


5.10 The transfer function of a system given by 

C(s) _ 10(s + l)(s + 3-)(s + 5) 

R(s) s(s + 2)(s + 4)(s + 5) 

indicates that this system is either uncontrollable or unobservable. Check by 
the rank test, after a zero-pole cancellation takes place, that the remaining 
system is both controllable and observable. 

5.11 A discrete model of a steam power system was considered in Mahmoud 
(1982) and Gajic and Shen (1993). 

(a) Using MATLAB, examine the controllability and observability of this 
system, represented by 


' 0.915 

0.051 

0.038 

0.015 

0.038 ' 

-0.030 

0.889 

-0.001 

0.046 

0.111 

-0.006 

0.648 

0.247 

0.014 

0.048 

-0.715 

-0.022 

-0.021 

0.240 

-0.024 

-0.148 

-0.003 

-0.004 

0.090 

0.026 


= [0.010 0.122 0.036 0.562 0.115 ] T 


C d 


110 0 0 
0 0 111 


(b) Find the system transfer function and justify the answer obtained in (a). 

5.12 Using MATLAB, examine the controllability of the magnetic tape control 
system considered in Chow and Kokotovic (1976) 


'0 

0.40 

0.000 

0.00 ' 


'0 ' 

0 

0.00 

0.349 

0.00 

, B = 

0 

0 

-5.24 

-4.65 

2.62 

0 

.0 

0.00 

0.000 

-10.00. 


.10. 
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5.13 Find the controllability matrix of the system in the phase variable canonical 
form and show that its rank is always equal to n. 

5.14 Linearize the given system at the nominal point (*i„, x- 2n , u n ) — (1,0, k) 
and examine system controllability and observability in terms of k 

x\ — x\u , aq(0) = 1 

x 2 — x\x 2 + u, «|(0) = 0 
y = x i 


5.15 Find the state space form of a system given by 

d 3 y d 2 y , dy du 

~dfi + W + k Tt Jr ' 2 y ^Tt +u 


and examine system controllability and observability in terms of k. Do they 
depend on the choice of the state space form? 

5.16 Given a linear system described by 


d 2 y 
dt 2 



+ y = 



u, 


2/(0) — o, 


dy{ o) 

dt 


Transfer this differential equation into a state space form and determine the 
initial conditions for the state space variables. Can you solve this problem 
by using an unobservable state space form? Justify your answer. 

5.17 Check that the matrix A given in (4.36) and the matrix C = \/Q defined 
in (4.39) form an observable pair. 

5.18 Derive formulas (5.55) and (5.56) for the reduced-order observer design. 

5.19 Using MATLAB, examine the controllability of a fifth-order distillation 


column considered in 

Petkov et al. 

(1986) 




'-0.194 

0.0628 

0 

0 

0 


1.306 

-2.132 

0.9807 

0 

0 

A = 

0 

1.595 

-3.149 

1.547 

0 


0 

0.0355 

2.632 

-4.257 

1.855 


0 

0.00227 

0 

0.1636 

-0.1625 

n 

[0 0.0632 0.0838 0.1004 

0.0063 1 T 


0 0 -0.1.396 -0.2060 -0.0128 
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5.20 Examine both the controllability and observability of the robotic manipula¬ 
tor acrobot whose state space matrices are given in Problem 3.2. 

5.21 Repeat problem 5.20 for the industrial reactor defined in Problem 3.26. 

5.22 Consider the state space model of the flexible beam given in Example 3.2. 
Find the system transfer function and determine its poles and zeros. Use 
Theorem 5.7 to check the controllability and observability of this linear 
control system. 

5.23 The system matrix for a linearized model of the inverted pendulum studied 
in Section 1.6 is given in Section 4.2.3. Using the same data as in Section 
4.2.3, the input matrix is obtained as 

B = [0 1 0 —2 ] T 

Examine the controllability of this inverted pendulum. 

5.24 A system matrix of a discrete-time model of an underwater vehicle is given 
in Problem 4.22. Its input matrix is given by Longhi and Zulli (1995) 

' 0.0258 -0.0002 0 ' 

0.102.3 -0.0010 0 

-0.0001 0.0258 0 

-0.0008 0.1023 0 

0 0 0.0055 

0 0 0 . 0221 . 

Check the controllability of this system. 



